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ABSTRACT 


This report presents the results of a study to investigate, by 
means of a computer simulation, the performance sensitivity of heli- 
copter IMC DSAL operations as a function of navigation system parameters. 
A mathematical model representing genetically a navigation system is 
formulated. The scenario simulated consists of a straight in helicopter 
approach to landing along a 6° glideslope. The deceleration magnitude 
chosen is ,03g. The navigation model parameters are varied and the 
statistics of the total system errors (TSE) computed. These statistics 
are used to determine the critical navigation system parameters that 
affect the performance of the closed-loop navigation, guidance and 
control systems of a UHIH helicopter. 
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INTRODUCTION 


Background and Objectives 

A truly unique characteristic of the helicopter is its ability to 
perform Decelerating Steep Approach and Landing (DSAL) maneuvers into a 
confined landing zone. In the past, such maneuvers have been performed 
only under Visual Meteorological Conditions (VMC) . That is , the helicopter 
pilot could perform a decelerated steep approach to hover only when he 
could see the target area from a sufficient distance to safely decelerate 
to a stop. Under Instrximent Meteorological Conditions (IMC) , the pilot/ 
autopilot must maneuver the aircraft along a desired flight path and 
deceleration profile so as to arrive at the Decision Height (DH) with a 
range rate that would permit a visual approach to a hover over the 
prescribed landing pad. 

A block diagram of a typical automatic. DSAL system is shown in 
Figure 1. The overall closed-loop system is relatively complex and con- 
sists of three distinct modules for performing the aircraft navigation, 
guidance and control tasks. The navigation system must process noisy 
position data (aircraft range, azimuth and elevation) provided by a 
ground-based landing guidance system (LGS) , using an on-board digital 
processor, to obtain estimates of the aircraft's position and range rate. 
These estimates are then used by the guidance and control systems to 
fly the helicopter along a desired position and velocity profile. 

The performance of the DSAL system depends upon a nxmiber of factors 
such as (1) commanded flight path and velocity profile, (2) aircraft 
dynamics, (3) on-board avionics and flight control sophistication, (4) 
on-board navigation filter characteristics, and (5) accuracy of the 
aircraft raw position information provided by the ground-based LGS. An 
investigation of all these five factors was considered beyond the scope 
of this effort. Instead this study is limited to an analysis of the 
DSAL system performance with varying on-board and ground-based (LGS) 
navigation system parameters. 
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The objectives of this effort are basically twofold: 

1. Conduct a generic sensitivity analysis of the helicopter DSAL 
system performance as a function of the on-board and ground- 
based navigation system parameters; and 

2. Investigate flight path and range rate total system errors 
during helicopter IMG DSAL tasks for various weather minima 
(DH « 200,150,100,50,0 ft). 


Technical Approach 

Due to the complexity of the DSAL system, this study uses a computer 
simulation to meet its goals. Specifically, an off-line computer program 
was developed to simulate the closed-loop DSAL system described by the 
functional block diagram of Figure 1. A nonlinear mathematical model is 
used to represent the UHIH helicopter dynamics [1] . Monte Carlo pro- 
cedures are used to obtain the ensemble means and standard deviations of 
several DSAL performance measures. The specific performance measures 
used in this study are described in the section on Simulation and Analysis. 
The details of the navigation and guidance laws and of the autopilot and 
flight control systems are given in the section on the Helicopter DSAL 
System and Appendix A, respectively. The on-board navigation filters are 
restricted to the commonly used a-0 filter structure [2] . The principal 
navigation signals affecting the helicopter's DSAL performance are range 
rate and range from the helipad center. Consequently, the scope of this 
study is limited to an investigation of helicopter DSAL performance using 
only range and range rate elements of the navigation system. Aircraft 
azimuth and elevation signals provided by the ground-based system are 
assumed to be perfect or without error. 
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HELICOPTER DSAL SYSTEM 


This section describes the specific DSAL guidance profiles and the 
mathematical models for the on-board (a-3 filter) and ground-based (LGS) 
navigation systems used in the off-line DSAL simulation program. 

The IMG DSAL maneuver involves decelerating to zero range rate while 
tracking the localizer and glideslope. Two deceleration guidance schemes 
have been investigated by several researchers: (1) Constant Deceleration 
Profile (CDP) [3] and (2) Constant Attitude Deceleration Profile (CADP) 

[4], This study investigated the performance of the navigation systems 
using CDP guidance and a 6° glideslope trajectory. 

DSAL Guidance Profile 

Figures 2(a) and 2(b) show the desired range rate and altitude profiles 
used in the simulation of the straight-in DSAL task. The nominal DSAL 
maneuver described by these profiles can be segmented into the following 
four distinct phases: (1) constant altitude (1500 ft) and constant air- 
speed (60 kts) cruise, (2) glideslope capture at 14271.5 ft followed by 
constant airspeed (60 kts) glideslope (6° G/S) and localizer tracking, 

(3) switch to constant range rate (60 kts) glideslope tracking at 8300 ft 

(4) glideslope tracking to a landing following initiation of a constant 
0,03g deceleration maneuver at 5285 ft. 

The constant deceleration profile provides a range rate guidance 
command (i^^) as a nonlinear function of range (r) defined by: 

r = / 2ar 

c 

where a is the deceleration magnitude. Figure 2(a) shows that after inter- 
cepting the deceleration profile the pilot has about 5285 feet to come to 
a touchdown at zero range rate. This deceleration range remains constant 
regardless of a glideslope angle at which the pilot approaches the helipad. 
However, the order in which the glideslope is captured and deceleration 
initiated depends on the approach altitude, selected glideslope, initial 



Conmanded Range Rate r (ft/sec) 




velocity and deceleration magnitude. In either case, the navigation, 
guidance and control systems enable the helicopter to fly along the glide- 
slope, localizer and parabolic velocity profile to come in to touchdown 
at zero range rate. 


Navigation System 

The azimuth, elevation, range and range rate signals that drive the 
helicopter flight control system and autopilot (described in Appendix A) 
are obtained by on-board digital processing of the raw position signals 
provided by the ground-based navigation system (LGS). In order to 
investigate the sensitivity of the helicopter DSAL performance to the 
navigation system characteristics, a mathematical model representing the 
acquired raw navigation signals from the ground-based LGS and the on- 
board processing of these signals is formulated as shown in Figure 3. 

The remainder of this section describes the details of the proposed 
mathematical model of the navigation system. Although, the model is 
general enough to describe the acquisition and processing of the range, 
azimuth and elevation signals, it was used only to simulate the range 
element of the navigation system. 

Groiind-Based Elements ; The ground-based LGS provides noisy measure- 
ments of range, azimuth and elevation at regular intervals. The noise in 
these measurements is modeled as a correlated random process with correla- 
tion time constant x seconds and standard deviation O . A constant bias 

n n 

B in the noisy measurement is also included. However, the effect of this 
bias on the DSAL performance is deterministic and can be removed through 
calibration. Hence, it’s effects are not considered in this study. 

Airborne Elements; In the on-board digital system, the raw measure- 
ment (y^) is sampled and quantized. The sampler has a frequency f and 
the quantization process Introduces a truncation (or roundoff) error de- 
fined by q units per Least Significant Bit (LSB) . q is a fimction of 
the full scale capability of the measurement and the word size of the 
digital processor. This discrete measurement is then filtered by the air- 
borne digital filter as shovm in Figure 3. 
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Figure 3. Navigation System Model 


The digital filter processes the raw (noisy) range measurements and obtains 
estimates of the helicopter's current range and range rate. Since the 
digital filter does not have a floating point processor, the estimated 
range rate is quantized as well. 


1 

j 
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The analysis below shows that the quantization of the raw range 

measurement (y ) contribute a mean and standard deviation to the discrete 
■^n 

time measurements y . Considering Figure 3 the raw signal y^ can be 
written as: 


n 


= y^ + N + B 


(1) 


n 
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r 


where N and B are additive colored noise and bias respectively. Further, 
this raw signal is quantized to give the signal [5]; 


n 


= y + e 


( 2 ) 


where y is the truncated signal and 

q 

e is the truncation error. 


Consequently, if the quantized signal y^ is measured then the residue r 
can be written as; 


r 


r = - yq 

= e - (N + B) 


(3) 




Assuming that the truncation errors are random with a uniform density [5] , 
and truncation and noise are independent processes, taking expected values 
of the above gives; 


r = - B 


2 q 2 

T = + a 

r 12 n 


where q is the quantization level of the range, 

r, are the mean and standard deviation respectively of the 
residue r. 


(4) 

(5) 


Range and Range Rate Filters; A critical part of the airborne 
navigation subsystem is the digital range/range rate filter. For the 
purpose of this study, a commonly used a-6 filter [2] was simulated. 
Appendix B describes the a-3 filter as implemented in this study. Figure 
B.l shows the configuration of this filter. The a-3 parameters define 
the bandwidth (co ) and damping (5) of the second order filter: 


n 


0 ) = / 3 

n 


( 6 ) 


g+g AT 
2 / 1 “ 


( 7 ) 


where AT is the dynamics update rate. 
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Using the mathematical models for the navigation, guidance and control i 

systems, a simulation was developed to study the sensitivity of the closed 
loop performance of the UHIH helicopter to the navigation system parameters. ~! 

The next section presents and discusses the results of this simulation. 
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SIMULATION AND ANALYSIS 


Modules were developed to simulate the navigation, guidance and 
control system models described above. These modules were combined with 
an existing program that simulated the nonlinear aerodynamics of the UHIH 
helicopter. The complete simulation program was then used to study the 
closed loop behavior of the helicopter DSAL system. 

The DSAL simulation begins following localizes capture at a range 
of 15000 ft. A straight-in approach simulation run consists of the 
following four phases: (1) cruise at constant altitude of 1500 ft and 
constant speed of 60 kts, C2) constant speed glideslope capture, (3) 
constant speed 6“ glideslope tracking until initiation of deceleration 
at 5815 ft range- to^go and (4) constant deceleration (.03 g) to land 
at the helipad at zero range rate. Initial helicopter cruise altitude 
is selected to be 1500 ft in order to minimize the effects of glideslope 
capture transients on the system performance measures at or near touch- 
down. The simulation run terminates when the helicopter altitude 
reaches ground level (0 ft). All simulation runs were done on the 
CDC 7600. 

System Performance Measures 

In order to evaluate the sensitivity of the overall system to the 
navigation system parameters, the performance measures shown in Table 1 
are used. Definitions (1-4) and (7-9) are measures of the Total System 
Error (TSE) , while definitions (5-6) are measures of the Navigation 
System Error (NSE) . 

NOTE: 

(a) Desired azimuth, elevation and altitude are computed in terms 
of the helicopter's actual range and desired glideslope. 

(b) Desired range rate is computed as a function of the helicopter's 
actual range. 

(c) Desired altitude at decision range is the altitude the helicopter 
should be at actual decision range on a desired glideslope path. 
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Table 1 


System Performance Measures 


No. 

Performance 

Measure 

. Definition 

Plot 

Variable 


1 

Localizer 

Deviation 

Desired Azimuth - 
Actual Azimuth 

LOG ERROR (DEG) 

— 

H 

Glides lope 
Deviation 

Desired Elevation - 
Actual Elevation 

GS ERROR (DEG) 


1 

Altitude Error 

Desired Altitude - 
Actual Altitude 

ALT ERROR (FT) 

"1 

■ 

Range Rate 
Error 

Desired Range Rate - 
Actual Range Rate 

RDOTD - RDOTA (FT/S) 

i 

j 

5 

Range NSE 

Actual Range - 
Estimated Range 

RANGE NSE (FT) 

i 

! 

6 

Ranga Rate 
NSE 

Actual Range Rate - 
Estimated Range Rate 

RDOT NSE (FT/S) 

i 

7 

Altitude Error 
at Decision 
Range 

(Desired Altitude - 
Actual Altitude) at 
Decision Range. 

Values of 2000, 1500, 
and 500 ft are chosen 
corresponding to Deci- 
sion Heights of 200, 
150, 100 and 50 ft, 
respectively on a 6° 
approach 

DECISION ALTITUDE 
ERRORS (y,a) FT 

“1 

1 

. 1 
”1 
n 

8 

Range Rate 
Error at 
Decision 
Range 

(Desired Range Rate - 
Actual Range Rate) at 
Decision Range. Deci- 
sion Range value as 
defined in 7 above. 

RANGE RATE ERRORS 
(y,a) KNOTS 

i 

9 

Range Error 
at Touchdown 

Desired Range at 0 ft 
Decision Height - 
Actual Range at 0 ft 
Decision Height. 

DECISION RANGE 
ERRORS (y,a) FT 
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Monte Carlo procedures with thirty individual runs are used to com- 
pute the ensemble mean and standard deviation of the performance measures 
defined above, The measurement noise sequences for each run are different. 
The ensemble statistics of variables (1-6) defined above are polotted as 
a function of range-to-go, The mean and one standard deviation envelopes 
(p + a) of variables 7 and 8 are plotted as a function of the decision 
range. 

i 

The remainder of this section is organized in the following four 
parts; 

1. Case descriptions of navigation system parameters, 

2. Overall system performance results for the DSAL system with 
(i) perfect range /range rate information and (ii) with the 
nominal navigation system. 

3. Sensitivity analysis results as a function of on-board 
navigation system parameters. 

4. Sensitivity analysis results as a function of ground-based 
navigation system parameters. 


Case Descriptions of Navigation System Parameters 


The navigation system mathematical model described earlier intro- 
duced eight parameters that characterize a navigation system. These 
variables are defined here again for convenience. Nominal values of 
these parameters for the range system are given in parentheses. 

LGS Parameters: 

= measurement noise : (1 ft) 

= noise correlation time constant: (0. 1 s) 

B = measurement bias:(0 ft) 

Airborne System Parameters: 

f = data update frequency = hz) 

q = quantization: (1 ft) 
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= rate quantization: (1 kt = 1.7 ft/s) 

C = a-3 filter damping: (0. 707) 

03^ = a-3 filter bandwidth: (2 rad/s) 

In order to investigate the effects of the navigation system para- 
meters on the overall system performance, a subset of combinations of the 
eight parameters defined above is selected. The study is conducted by 
varying some of the parameters defined above about a nominal set. 
Variations of B and C are not considered. Any effects due to the 
changes in can be minimized by an appropriate selection of the filter 

bandwidth O) . Effects due. to the bias B are deterministic and can be 
n 

corrected for. The filter damping ? was fixed at .707. Table 2 below 
shows the combinations of the navigation system parameters investigated 
in this study. 


Table 2. Navigation System Parameters Used in Study 


Case 

On-board System Parameters 

Ground-based 

System 

Parameter 

a„(£t) 

f(Hz) 

q(ft) 

q^(ft/s) 

03 (rad/s) 
n 

0 (Nominal ) 

16.0 

1.0 

1.7 

2.0 

1.0 

1 

4.0 

- 

- 

- 

- 

2 

8.0 

- 

- 

- 

- 

3 

- 

10.0 

- 

- 

- 

4 

- 

30.0 

- 

- 

- 

5 

- 

- 

8.5 

- 

- 

6 

- 

- 

17.0 

- 

- 

7 

- 

- 

- 

0.2 

• 

8 

- 

- 

- 

10,0 

- 

9 

- 

- 

- 

- 

10.0 

10 


- 

- 

- 

30.0 


Note that entries with a dash (-) mean that the nominal value (case 0) 
was used for the pertinent parameter. 
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Cases 1 through 8 are used to investigate the performance sensitivity 
to on-board navigation system parameters (f,qjqy,w^). Case 9 through 10 
show effects of varying the ground-based LGS parameter (<^^) • 

The next section compares results of the overall system performance 
using perfect range/range rate information and the notainal navigation 
system. These results will then be used as a reference to evaluate the 
sensitivity of the overall closed-loop system to individual navigation 
model parameters. 

System Performance Using Perfect Range/ 

Range Rate and the Nominal Navigation System 

In order to separate tracking errors caused by the control and 
navigation systems, the simulation was first run with perfect range and 
range rate information (i.e,, a perfect navigation system is assumed 
with zero range measurement noise and no on-board quantization or filtering) . 
It was found that when the range rate command as shown in Figure 2(a) is 
used the helicopter touches down short of the helipad with a non-zero 
range rate. In order to decrease this touchdown velocity, a 5 ft/s 
velocity offset was subtracted from the nominal commanded velocity of 
/I, 96r. This velocity offset meant that deceleration began at 5815 ft 
instead of the nominal 5285 ft. 

Figures 4(a)-(h) show the ensemble statistics of performance measures 
(1-6) of Table 1 using (i) perfect range/range rate and (ii) the nominal 
navigation system parameters (Case 0 in Table 2). Figures 4(a)-(c) show 
that the localizer, glideslope and altitude deviations are negligible. 

These errors are large at close range due to the ill-defined nature of 
these variables at the origin of the reference frame. Figure 4(d) shows 
that both simulations terminated with the helicopter landing short at a 
range of 28 ft and moving at about 3 ft/s (1,75 knots). This was re- 
garded as an acceptable touchdown. Since touchdown is defined to occur 
when the helicopter reaches ground level, compensation must be made in 
the vertical guidance loop in addition to the range rate loop to ensure 
that the simulation terminates with the helicopter at zero range and 
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Figure 4. Ensemble Statistics for DSAL System with Perfect Range/ 
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Figure 4. Ensemble Statistics for DSAL System with Perfect Range/ 
Range Rate and Nominal Navigation System 
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range rate. Although compensation for range rate was made, the vertical 
guidance was not modified. Hence results with the nominal navigation 
system are used as a baseline to study performance sensitivity to naviga- 
tion system parameters. 

Figure 4(e) shows about .5 ft mean range NSE, This corresponds to 
the expected mean range NSE due to the range quantization (q/2 in Eq. (4)). 
Similarly, Figure 4(f) shows a mean range rate NSE of 0,85 ft/s due to 
range rate quantization (q^/2 in Eq, (4)). 

The altitude error and range rate error at decision range are not 
shown, since these errors are very small. 

Sensitivity to On-board Navigation System Parameters 


The on-board receiver and range/range rate filter are characterized 
by the following four parameters: 

(a) f(Hz) - range measurement sampling frequency, 

(b) q(ft) - range measurement quantization level. 

(^) q,j.(ft7s)- quantization level of the estimated range rate. 

(d) CO (rad/ s) - bandwidth of second order (a-8) range/range rate 
filter. 

The sensitivity of the overall system performance to the above para- 
meters was studied by varying them one at a time while holding all others 
at the nominal value (see Table 2) . The essential results are presented 
in Figures 5(a)-(c) which show plots of decision altitude, range rate at 
decision range and of range TSE at touchdown, respectively. The following 
observations are pertinent; 

1. The closed-loop system performance is not very sensitive to 
range measurement sampling frequency (f) and quantization (q) . 

2. Increasing the quantization of the estimated range rate causes 
oscillations in the system as shown by the decision altitude 
errors (Figure 5(a)). Coarser range rate quantization 
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Increases the magnitude of the mean and standard deviation 
of the range rate TSE (Figure 5(b)). The standard deviation 
of the range rate TSE at touchdown is about 7 knots when a 
range rate quantization of 10 knots (17 ft/s) is used. Con- 
sequently, the helicopter flight profile is higher and faster 

than desired and displays large run-to-run variability during 
the final 500 ft prior to touchdown. 

3. Decreasing the bandwidth of the range/range rate filter to 

,2 rad/a causes oscillations in the system as shown by Figures 
5(a),(b), These oscillations are more evident in the ensemble 
statistics of the TSE as a function of range-to-go (Figures 
9(c). Cd)). 

Figures 6(a)^(h) show the ensemble statistics as the range sampling 
frequency (f) is varied. These figures confirm that the DSAL perfor- 
mance is not very sensitive to f. 

Figures 5 (a) , (b) show that the altitude and range rate TSE are 
insensitive to range quantization (q). However, these plots convey no 
information at ranges less than 500 ft. Figures 7(a)-(h) show the 
performance sensitivity to range quentization (q) . Figure 7(d) shows 
that the range rate TSE increases as the helicopter comes in to touch- 
down. Note that the performance of the nominal system (solid line in 
figures) is . to be used as a basis for comparison. Figure 7(c) shows 
that the effects on altitude TSE as q is increased are not noticeable. 
Figure 7(e) shows that the range quantization causes a bias in the range 
NSE as shown by equation (4) . It is this bias in estimated range that 
causes a corresponding error in the range rate TSE. 

The oscillations in the helicopter flight profile due to range rate 
quantization are evident in Figures 8(a)-(h). Once again, note the bias 
introduced in the range rate NSE (Figure 8(g)) due to the range rate 
quantization. The figures also show that the mean and standard deviation 
of range rate NSE decrease drastically around 8000 ft and then increase 
again. This step change occurs when the pitch axis control system 
switches to range rate feedback Instead of airspeed at about 8300 ft 
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Figure 5(b). Sensitivity of Range Rate TSE to On-board Navigation Parameters 
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range-to-go. Since the aircraft is commanded to maintain a range rate 
of 101,0 ft/s, the range rate quantization of 8.5 or 17.0 ft/s would 
cause the estimated range rate to be 93,5 or 85 ft/s, respectively. 
Consequently, this results in the large range rate NSE for ranges 
greater than 8300 ft. However, when the range rate is fedback instead 
of airspeed, the aircraft speeds up and so its actual range rate is 
slightly greater than 102,0 ft/s, With this speed, quantizations of 
8.5 or 17 ft/s result in the estimated range rate being 102.0 ft/s 
and therefore close to the true range rate. Consequently, the range 
rate NSE decreases until the aircraft begins decelerating at about 
6000 ft. 

The digital a-8 filter has to be designed so that (i) the range 
measurement noise suppression characteristics are good and (ii) lags 
introduced in the filtered data are not so severe as to cause system 
instabilities. Such a design calls for a compromise in the filter band- 
width (6J^) . Sensitivity of the DSAL system to this parameter (co^) can be 
seen in Figures 9(a)-(h). Decreasing the filter bandwidth (co^) from 2 rad/s 
to 0.2 rad/s causes an oscillation in the system as shown by Figures 9(c) , 
(d) . The standard deviation of the range NSE decreases as shown by Figures 
9(f) as expected because of the lower filter bandwidth. If the bandwidth 
is increased to 10 rad/s, the mean TSE are marginally different from the 
nominal system. Figure 9(f) shows that with the increased bandwidth of 
the a-8 filter, more of the range measurement noise "passes" through the 
system. However, since this noise is very small (a^ = 1.0 ft) , the effects 
on the TSE are minimal. 

The results of the simulation have shown that the DSAL system is 
most sensitive to estimated range rate quantization and filter bandwidth. 

An increase in the former and a decrease in the latter cause system oscilla- 
tions. A decrease in range measurement sampling frequency (f) to 4 Hz 
causes no degradation of the system performance. Range quantization effects 
are felt at short ranges-to-go. An increase in this quantization results 
in a larger range rate TSE. 
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Sensitivity to Ground-Based Navigation 
(LGS) System Parameter 


The ground-based segment of the navigation system has been modeled 
such that the range measurements are assumed to be corrupted by colored 
noise. The single parameter used in this study to characterise the 
noise properties of the range measurements is the standard deviation a 

n 

of the noise process. 

Figures 10(a)- (c) show the effects of increasing the range noise 

a from 1 to 30 ft. Note that a three fold increase in the range noise 
n 

standard deviation from 10 to 30 ft, causes almost a ten fold increase 
in the standard deviation (a) of the decision altitude TSE. Figures 
10(a)-(c) show that increasing the range noise standard deviation 
causes an increase (in magnitude) of the total system errors. This 
effect is due to the non-linear nature of the overall system being 
simulated. 

The ensemble statistics as a function of range-to-go when the noise 
level (a^) is varied are shown in Figures ll(a)-(h). Figures 11(c) and 
(d) show that with a of 30 ft, the mean altitude and range rate TSE 

are poorer than with = 1 or 10 ft. Note that Figure 11(d) shows that 

deceleration begins earlier when the noise level increases. This can be 
explained by considering Figure 11(g). The decision to initiate decelera- 
tion is made by comparing the commanded range rate and the estimated 
range rate. The wide excursions in the mean range rate NSE imply that 
over certain segments of the flight the estimated range rate is larger 
than the true rcinge rate causing an earlier deceleration. 

These simulation results show that: (a) the increases in the stand- 
ard deviation of the range noise are felt in the mean of the TSE due to 
the non-linear nature of the system, and (b) the system performance de- 
grades very rapidly as the range noise level is increased. 
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The standard deviation of the range noise is one measure of the 
accuracy of the LGS and reflects on the price/performance tradeoff of the 
ground-based system. Since the impact of the range noise on the DSAL 
system performance is significant, it is important to study the system 
sensitivities to the on-board navigation parameters using a larger nominal 
range noise (a^) . 

Cases 1-8 of Table 2 were repeated using a nominal of 10 ft. 

Figures 12(a)-(c) show the system performance sensitivity (measures 

7-9, Table 1) as a function of on-board navigation parameters with = 10 

ft. Comparing Figures 5(a) and 12(a), it can be seen that the system 

performance is not very sensitive to range measurement sampling frequency 

for both a “ 1 and 10 ft. 
n 

Although range quantization does not seem to have an effect on the 
decision altitude or range rate TSE (Figures 5(a),(b)) with a nominal 

- 1 ft, the effects on these performance measures are noticeable with 

- 10 ft and q = 30 ft (Figures 12 (a), (b)). Note that the effects of 
the range measurement noise and range quantization on decision altitude 
TSE are not merely additive (Figure 12(a)). Further, the standard devia- 
tion of the decision altitude errors increases as the helicopter comes in to 
touchdown when the range quantization level is made larger (Figure 12(a)). 

The effects of range rate quantization when = 1 or 10 ft are 
similar; both systems show oscillations when a range rate quantization 
(q^) of 17 ft/s is used. 

Figures 5(a), (b) and 12 (a), (b) both show that decreasing the band- 
width of the a-B filter to .2 rad/s causes the helicopter system to 
oscillate. On the otherhand, increasing the bandwidth of the range/range 

rate filter to 10 rad/s when a = 10 ft did not allow the DSAL simulation 

n 

to terminate normally. This is because the larger filter bandwidth causes 
the range noise to "pass" through resulting in unacceptable range and 
range rate estimates. 
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A comparison of Figures 5(c) and 12(c) indicates that increasing the 
range measurement noise increases the magnitude of the mean range touch- 
down error. This was observed in Figure 10(c) as well. Figure 5(c) shows 
that as the range rate quantization is increased, the magnitude of the 
mean touchdown error increases. However, if a a of 10 ft is used, then 
this error decreases as the range rate quantization is increased (Figure 
12(c)). This seems to indicate that the variations in the touchdown 
range TSE due to range measurement noise and range rate quantization are 
in opposite directions. This is not to say, that one can have a very 
large range measurement noise and compensate it by a very large range 
rate quantization. This is because the helicopter DSAL system oscillates 
a's the range rate quantization increases. 
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CONCLUDING REMARKS 


A closed loop computer simulation of the UHIH helicopter DSAL system 
was developed for investigating performance sensitivity to navigation 
model parameters. A navigation model consisting of on-board and ground- 
based (LGS) elements with five (f parameters was formulated. 
The first four parameters characterize the on-board system and the last 
one characterizes the LGS. These parameters were varied about a nominal 
set to study the sensitivity of the helicopter DSAL performance. Perfor- 
mance measures were defined in terms of the total system errors of the 
following flight path variables: localizer, glideslope, altitude, range 
rate, altitude and range rate at decision ranges and range at touchdown. 


The following conclusions on system performance parameter sensitivity 
were drawn from simulation results: 

1. On-Board Sampling Frequency (16, 8, 4 Hz). 

System performance was not degraded by decreased sampling 
frequency. 

2. On-Board Range Quantization (1, 10, 30 ft). 

System perfonnance was not degraded for increased range 
quantization, except only slightly when increased to 30 ft 
for the nominal case with ground-based range noise equal 
to 10 ft. 

3. On-Board Range Rate Quantization (1, 5, 10 kts) . 

Increased range rate quantization caused an almost propor- 
tional increase in negative bias for range rate error. 

Increased range rate quantization also caused system 
oscillations to increase. 

4. On-Board Filter Bandwidth (0.2, 2.0, 10.0 rad/sec). 

A filter bandwidth of 0.2 rad/sec was too low and caused an 
undesirable increase in glideslope tracking oscillations. A 
filter bandwidth of 10.0 rad/ sec caused the system to go 
unstable when ground-based range noise was increased from 1 
to 10 feet, 

5. Ground-Based Range Noise (1, 10, 30 ft). 

System performance was significantly degraded when ground-based 
range noise was increased from 10 to 30 feet. When ground-based 
range noise was increased to 30 feet , along-track range error at 
touchdown (-70 + 65 ft) was significantly greater than in all 
other parameter sensitivity cases tested. 
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APPENDIX A 


DSAL AUTOPILOT/FLIGHT-CONTROL SYSTEM 

This appendix describes the four-axis stability augmentation system 
(SAS) and autopilot used in the simulation of the overall closed-loop DSAL 
system. 


Collective Axis 

The autopilot flight path coupler in the collective axis operates in 
the altitude hold mode (ALT) or in the approach mode (APP) (Figure A.l). 

Altitude hold is accomplished by synchronizing to the aircraft's 

altitude at the time of engagement and then applying the error signal as 

a command to the collective axis. The error signal is augmented by the 

error rate which is obtained by dif ferentation (with a lag) of the error. 

This error rate is complemented by attitude (9) stabilized vertical 

acceleration (a ) . 

z 

In the approach mode (APP) , the filtered range (r) and angular glide- 
slope deviation are multiplied to give vertical deviation about the desired 
glideslope. This error signal is augmented by a computed error rate (as 
in the altitude hold mode) . The error plus the error rate signals then go 
through an integrator to provide a proportional plus integral plus deriva- 
tive (PID) controller. The collective stick position is fed forward 
through a 1.5 second washout filter to the PID signal. The above summed 
signal is then applied to the collective axis autopilot (Figure A. 2) . 

The collective axis autopilot has no SAS and therefore the flight path 
coupler command goes straight through the series and parallel servo to 
generate the collective axis command. 
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Figure A. 2. Autopilot Collective Axis 
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Pitch Axis 


The pitch axis uses a relatively simple control law to control air- 
speed, or range rate (Figures A, 3, A, 4), It operates in one of two modes; 
speed hold (SPD) or deceleration (DECL) , In the speed hold mode, the 
synchronizing speed is compared with the airspeed. The airspeed being 
rather noisy is filtered by a .1 second lag first order filter. The 
filtered airspeed is complemented by attitude (9) stabilized longitudinal 
acceleration (a ). The error signal is then limited and added to a washed 
out collective stick position signal. The total signal is then multiplied 
by a gain of .57 deg/ft/sec to give the pitch command (6^)» At a range 
of 8300 feet or less, the DECL mode goes from ARMED to ON and the esti- 
mated range rate is fed back instead of airspeed. However, the actual 
deceleration does not begin until the deceleration profile is intercepted 
(approximately at 5285 feet for a cruise speed of 60 knots and a decelera- 
tion of .03G). The deceleration profile generator computes the commanded 
range rate as a function of estimated range. (A constant velocity offset 
is subtracted from the commanded range rate to reduce the range rate error 
noticed in the simulation). The commanded range rate is continuously 
compared to the aircraft's speed. In the speed hold mode the resulting 
etror signal is limited to zero so that no speed up command is generated. 
When the deceleration profile has been intercepted (DECLINT) , the error 
signal is limited to + 5 ft/sec. As in the speed hold mode, the error 
signal is added to washed out collective stick position signal. The total 
signal then goes through a gain of .57 deg/ft/sec to give the pitch command 
9^ which then acts as an input to the pitch axis SAS (Figure A. 4). 

The pitch axis SAS is veiry similar to the roll axis , described in the 
following section. 


Roll Axis 

Figures A. 5, A. 6 show the roll control system. The roll axis commands 
the aircraft towards the localizer. The lateral displacement is computed 
by multiplying the filtered range (r) and angular localizer deviation (a^) . 
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The filtered range (r) is obtained by onboard processing of the DME signal 
from the ground based navigation system. This filtering of the range 
signal is described under Navigation Systems. The lateral displacement 
rate is derived from the computed displacement. The rate has a 5 second 
lag built in it because of the noise in the computed lateral displacement. 
This is complemented with attitude (cj)) stabilized lateral acceleration (a^) . 
The lateral displacement and the lateral rate (complemented as above) , are 
added to the integral of the localizer angle to form a proportional plus 
integral plus derivative (PID) type controller. The above signal is then 
multiplied by a gain and limited to form the roll command . The roll 
command is then fed into the stability augmentation system (Figure A. 6). 

Roll attitude and roll rate are fed back to produce a lateral cyclic command 
through a series and a parallel servo. The roll rate is obtained by 
differentiating the roll angle assuming no rate gyros are available. 


Yaw Axis 


Figure A. 7 shows the stability and control augmentation system (SCAS) 
for the yaw axis at high speeds (turn coordination mode). Since the air- 
craft tends to yaw when the collective is dropped or raised, the collec- 
tive stick position (radians) acts as the primary input to the system. 

The turn coordination mode requires lateral acceleration, yaw and roll 
rate feedbacks as shown. The pedal command is generated by passing the 
yaw signal through a series and parallel servo. In the heading hold 
mode, the commanded heading angle ( 4 >^) acts as an input in addition to the 
collective stick position. This mode requires only the yaw angle and 
yaw rate feedback (Figure A, 8). 
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Figure A. 7. SCAS: Turn Coordination Mode 















APPENDIX B 


RANGE AND RANGE RATE FILTERS 


g-B Range/Range Rate Filters 

The range/range rate filters investigated in this study is the 
classical second order g-B filter [2], Figure B.l shows the analog 
representation of this filter. In this filter r^ is the range measure- 
ment, and ^2 estimates of the range and range rate res- 

pectively, The state space representation of the filter is: 



Figure B..1, Analog Representation of the g-B Filter 
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Equations B..2 and B..3 show that the filter bandwidth and damping ratio 
are related by: 


Discretizing B.l by a simple Euler explicit scheme gives the following 
discrete equations for and X 2 : 


(n+1) 


$ 2 (n+1) 


1 - oAT AT Xj^(n) 


-bat 


X2(n) 


r (n+1) 
in 


where AT is the measurement update period. This scheme requires: 

|oAt1 << 1 E.7 

for the discrete solution to be accurate. Consequently, the measurement 
update period has to be sufficiently small to obtain accurate estimates of 
range and range rate. 

Since the Of-B filter is based on a kinematic model where the acceleration 
is assumed zero, the filter can be implemented in a different way in which 
the measurement update period need not be small. The kinematic model can be 
written as : 


0 1 


0 O x 


= A X 


r^ = [1 0] 

m 


= h X 


where x^^ and X 2 are the range and range rate of the helicopter. Using a 
Kalman filter representation for the model described by B.8, B.9 gives the 
following procedure for estimating x^^ and X 2 : 
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B.IO 


X (n) = x(n) + K(r (n) - x(n)) 

— — m — — 

X (n+ 1 ) = <I> x(n) 


$ 



AT 

1 


K 


"a . AT 
3 • AT 


B.ll 


B.12 


B..13 


B.IO is termed the measurement update equation and E.ll is the time update 
equation. In between measurements the range (x^) and range rate (X 2 ) esti- 
mates are propagated by use of equation B.ll and when a new measurement is 
available, the innovation is added to the propagated estimates as shown in 
equation B.IO, 


The steady-state performance of the filter represented by equations E.IO 
through B.13 is best studied in the Z transform domain. Equations B.IO 
through B.13 can be written as; 


x^(n+l) = Xj^(n) + AT * X 2 (n) + a . AT (r^(n+l) - x^(n) - AT • X 2 (n)) 

B..14 

X 2 (n+ 1 ) = ^ 2 (n) + 3 • AT (r^(n+l) - Xj^(n) - AT • X 2 (n)) 

E-15 


Taking Z transforms of B.14 and B..15 gives: 


where 


%(z) 

r (z) 
ni 



z r(z-l)g + 3 at] 

(z-1)^ + (a' + 3 ’ at) (z-1) + 3 ' at 

a • AT 

3 • AT 


B.16 


If the error e(t) is defined as; 
e(t) = r^(t) - 


then the transfer function e(z)/r (z) can be obtained from B,16; 

m 
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K.17 


e(z) 

t (z) 

in 


(z^l)^ g-q") 

A(Z) 


n 


where 


A(z) = Cz-1)^ + (a + b' at) (z-1) + b' at 


B.18 


Using the final value theorem 16], the following equations show the steady- 
state behavior of the error e(t): 


For a step r (t), Lim e(t) 
m 


For a, ramp 


Lim (z-1) e(z) 
z 1 


Lim (z-1) 
z 1 


(z-1)^ (1 
A(z) 


I 

-a ) 2 - n 

■ (z-1) “ 


Lim e(t) 


Lim (z- 
z ^ 1 


[(z-l)^ (l-a')l ATz 


R.19 


= 0 


B ..20 


For a parabolic r (t), Lim e(t) = Lim (z-1) . 

m * 


= ^ • AT 


B.21 


(z-1)^ d-a') 

(AT)^ z(z+l) 

n 

L A(z) J 

z(z-l)^ 

i 


For the DSAL law used in this study, the range is a parabolic function 
of time. Consequently the a-B filter will show a steady state error in 
the range estimate as shown by equation B.21. 


Equations B.14 and B.15 can also be used to derive the transfer function 
X 2 (z)/v(z) where v(z) is the Z transform of the velocity function v(t) : 


X2(z) 

v(z) 


B z . AT 
A(z) 


B.22 


e (z) 
v(z) 


(z-1) ((z-1) + g ) 
A(z) 


B.23 
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For a step velocity input the steady-state error can also be obtained 

by using the final value theorem; 


Lim e (t) 

V 


Lim (z-1) 
z-*-l 


(z-1) 


((z-1) + g ) 
A(z) 


z 

(z-1) 


= 0 


B..24 


For a ramp input velocity, the steady state error is: 


Lim e (t) 

V 


Lim (z-1) (z-1) 
z->-l 


((z-1) + g ) 
A(z) 


AT 

(z-l)‘ 


a' 


B.25 


The parameters a and 6 are related to damping and bandwidth of the 

analog filter via equations B.4 and B.5. However, since the filter is 

implemented digitally, the mapping of a and B into 5 and co has to be 

sAT 

derived usihg equation B.16. Since z = e , substituting B.26 for z in 
B,16 gives B.27 below. 


z = 1 + 

sAT 


B.26 

x^(s) 

(1+ sAT) (as + B) 


B.27 

t (s) 
m 

s^ + (a+ BAT)s + B 


Assuming that 

1 + sAT = 1, B.27 can be 

reduced to : 


x^(s) 

(as + B) 


B.28 

r (s) 
in 

s^ + (a + BAT) s + B 


Equation B.28 

shows that the bandwidth 

(co^) and damping (?) are 


related to a and B via B.29 and B.30: 



B = CO ^ 
n 



B.29 

a = 2 ? 0 )_ 

- bat 


B.30 


n 
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It can be seen that B.29 Is identical with B,4 whereas B,30 is equal 
to B.5 only if 

AT « ^ . 

CO 
n 


B.31 



, ) 


n 

i 

n 

i 



] 
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